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Abstract 

Properties of superfluid states of two-dimensional electron systems 
with critical antiferromagnetic fluctuations are investigated. These 
correlations are found to result in the emergence of rapidly varying 
in the momentum space terms in all components of the mass oper- 
ator, including the gap function A(p). It is shown that a domain, 
where these terms reside, shrinks with the temperature, leading to a 
significant difference between the temperature Tc, at which supercon- 
ductivity is terminated, and the temperature T*, where the gap in the 
single-particle spectrum vanishes. 

The problem of high-temperature superconductivity is known to defy so- 
lution within the Fermi liquid approach. Initially it looked like it had to 
do with only normal states of high- Tc superconductors. However, later it 
was acknowledged that the BCS theory fails in dealing with their superfluid 
states, as well [H This is best demonstrated by the persistence of a gap 
A in the spectra of single-particle (sp) excitations of many high-T^ super- 
conductors above the critical temperature Tc, at which superconductivity 
disappears (the so-called pseudogap phenomenon IH El El HI)- 

Another salient feature of two-dimensional electron liquid of high-Tc su- 
perconductors is the universality of its phase diagram versus the doping x. 
At low < Xc — 0.05, corresponding to the filling, close to 1/2, two- 
dimensional compounds are antiferromagnetic insulators. At larger x, an- 
iferromagnetic ordering is nil, but in the vicinity of the phase transition, 
long-range correlations with wave vectors q, close to the antiferromagnetic 
vector Q = (7r,7r), turn out to be drastically enhanced, which results in the 
divergence of the electron-electron scattering amplitude F = Fq + Ta(Jia2 
with 

F,(q ^ Q, ^ 0; a;) ~ [(q - Q)^ + r'^x) + ic\u;\]-' , (1) 
the correlation radius ra{x) becoming infinite at x = Xc ISHEl- 
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The impact of this singularity on sp properties is studied proceeding from 
the RPA formula = (Fq * G), which presents an associated with antifer- 
romagnetic fluctuations part of the mass operator S as a convolution 
of the amplitude Fq with the sp Green function G ( see e.g. [5-9]). For 
a long time, attention was focused on the energy dependence of S, while 
its momentum dependent part S(p, e = 0) was parameterized by the ef- 
fective mass m*. This is justified in systems with short-range correlations, 
where the mass operator S(p,£: = 0) is a smooth function of p. But this 
is not the case. Straightforward calculations show that long-range correla- 
tions trigger a rapidly varying with p component Sa(p) of the function 
E(p, £ = 0) = Sr(p) + Sa(p), being a convolution of F^ and the pole part 
G"' of the Green function G. It should be emphasized that Sa(p) has to be 
evaluated self-consistently, otherwise the flattening of the sp spectra ^(p) in 
normal states, found in Ref. JU] and observed in many high-Tc compounds, 
gets lost. 

To get rid of the energy-dependent terms in S we calculate the derivative 
dReTi{p,e)/dp (ReF^ * dlm.G'^/dp). After simple algebra we obtain 

= -zj Up -p,,uj = E{p))^^dn . (2) 

Here z = [1 — is the renormalization factor, dr = d'^p/{2'nY, and 

n(p,T) = v\p){l - f{E)) + (1 - v\p))f{E) = \ - ^tanh:^ (3) 
is the quasiparticle momentum distribution. In this formula, f{E) = (1 + 



eME/T)r\ while v\p) = {E{p) - ^(p)) /2E(p), where E(p) = ^^(p) + ^Hp) 
and A is the gap function, while ^{p) = z + E„(p) + ^^(p)) = ^°(p) + 
zJ]a{p) is the sp energy spectrum of the normal state measured from the 
chemical potential /i. To a good approximation, the spectrum ^°(p) and the 
LDA electron spectrum ^jj are related by ^°(p) = ^p/m*. 

In what follows, the argument u = E{pi) of the function Ta{p — Pi,^^^) 
in Eq. (j2I) is replaced by 0, since both the functions dn{pi,T)/dpi and 
^a{p — Pi), taken at ^(p) = 0, are peaked at ^(pi) = 0. Upon inserting this 
result into Eq.(j21) and integrating over momenta, one finds 

e(p) = e°(p) + J Taip - Pi)n(pi, T)dn . (4) 
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The gap A(p), obeying the BCS gap equation, is also decomposed into a 
sum A(p) = Aa(p) + Ar(p) of a regular Ar(p) and a rapidly varying with p 
component Aa(p). In the case of singlet pairing, the respective equation for 
Aa{p) reads (llj : 

r tanh ^^P^^ 

A,(p) = -3z' J r.(p - Pi)^^^A(Pi) dn . (5) 

The analysis of solutions of the above nonlinear equations is greatly facil- 
itated, if the interaction ((T)), taken at = 0, is approximated by a 5-function 
^z'^Ta{q) —>■ fa^i^l — Q) 05 appropriate in a domain of the momentum 
space, where the functions n(p) and Aa{p) change slower, than the ampli- 
tude Ta{p — pi). As a result, integrations cancel, and we are left with 

e(p) = e°(p) + fan{p - Q, T), e(p - Q) = e°(p - Q) + fanip, T) , (6) 

A. N ^xtanh(E(p-Q)/2T) ^, ^, ^ ^ , . tanh(E(p)/2T) 

A(p) = -/.A(p-Q)^^|^-^^, A(p-Q) = -/.A(p) 

(7) 

where the constant fa > 0, is small compared to the band width uq. In 
obtaining these equations we neglected the term Ar{p). Setting here T = 0, 
we arrive at a set of equations derived in Ref. ^2] in a different way. 

If pairing correlations are somehow suppressed, Eqs.Q are knocked out. 
Upon solving the two remaining Eqs.® we find that in case the van Hove 
points (±7r, 0) and (0, ivr) are situated quite close to the Fermi line (FL), a 
portion of the sp spectrum, adjacent to the van Hove points (vHP), turns 
out to be flat [12]. We shall see later that the flattening holds, if pairing 
correlations come into play. 

As seen from Eq.([7I), the gap A(p) changes its sign going over to a 
neighbor vHP, as in the conventional D-pairing model, in which the gap 
A.d{p) ~ {cospx — cospy) 7^ anywhere in the momentum space but the 
zone diagonals. However, in contrast to this model, nontrivial solutions of 
Eqs.(|7j) exist only in a domain C, boundaries of which are found by combining 
two Eqs.©, which yields 

,, taiih(i;(p)/2r) tarih(E(p - Q)/2T) 

— m Ei^) — 

Otherwise A^ = 0, and E{p) = |^(p)|, as in the Nozieres model [T3 | IT U IT5] 
with the effective long-range interaction r(q) ~ 5(q)- 
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In overdoped compounds, the domain C is made up of two quite narrow 
stripes. The first, denoted further Cp and described by equation ^°(p) = 0, 
is adjacent to the FL. The second, associated with the conjugate hne (CL), 
is determined by equation ^°(p — Q) = 0. In these compounds, the FL and 
CL are well separated, and when dealing with p G Cp the energy E{jp — Q) 
can be replaced by |^°(p — Q)|, so that Eq.® is recast to 

^(P) = 9{V) tanh(E(p)/2T), p G , (9) 

with the coupling constant g{j>) = fa/\^'^{p — Q)|- 

As X drops, the FL and the CL approach to meet each other at a critical 
doping Xm- In most of high- Tc compounds, such as Bi2212, Bi2201 etc., 
the FL is concave, while the CL, respectively, convex, and the first meeting 
between these lines occurs at the vHPs. Close to the vHPs, boundaries 
of the C domain are calculated combining Eq.(jSI) with relations E{p) = 
|^°(p)| and -E'(p) = |^°(p) + /a|, respectively, which yields restrictions —2/ < 
,^°(p) < /. In this case, Eq.® is easily solved, and close to the vHPs, the 
sp spectrum turns out to be quite flat: E{p,T = 0) ~ fa- We see that this 
value is significantly in excess of those obtained, if the FL and the CL have no 
points of intersection. When the gap landscape is drawn in the doping region 
X ~ Xm, it comprises four "twin towers", each one being associated with its 
own vHP. Each tower, whose height A^hp(T = 0,Xm), according to Eq.(jH)), 
equals /„, is connected with its neighbors by narrow "walls". According 
to Eq.(jni), their height drops towards the zone diagonals, where the gap A 
vanishes. Thus, we infer that the gap function A{px,py) attains its maximum 
AmaxiT = 0) = fa at the vHPs. This picture, confirmed by numerical 
calculations of Ref. ^21, is in agreement with the available experimental 
data 0. 

As T rises, the region C, where Aa(p) ^ shrinks, the effect, found 
first in Ref. ^3] under investigation of the Nozieres model JHl- Indeed, for 
points, fairly far away from the vHPs, Eq.ijU)) can be employed. Its nontrivial 
solutions exist only if g{p) > 2T. Since the function |^°(p — Q)|, identifying 
the energy splitting between the FL and the CL, rises, while the magnitude 
of the function g{p) drops, respectively, as the vector p moves along the 
FL towards the zone diagonal, the shrinkage begins in the diagonal region at 
Ti{x) ~ fa/ (2|^ma^(x)|, where ^max{x) is the bare sp energy, corresponding to 
the point of intersection between the CL and the zone diagonal. With further 
T increase, the shrinkage region is augmented, approaching the vHPs, where 
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the gap A has its maximum value. Eventually, the whole C domain shrinks 
into several symmetric points at the FL, closest to the vHPs. Recently, 
such a behavior of the gap landscape, the so-called arc phenomenon [2] was 
experimentally observed. The final shrinkage temperature T*{x ~ Xm) is 
easily evaluated from Eq.®. It is T*{x ~ Xm) — fa/'^, so that the gap 
^max{T = 0) and T* are connected with each other by 

A^aAT = 0)-2T* , Xr^Xm, (10) 

being in accord with the available experimental data P3I21- 

It is worth noting that inside the C region, the behavior of A remains 
the same as that in the BCS theory, since by retaining in Eq. ^ the leading 
terms, one obtains 

A2(p;T ^ T*) ^ 12T*(T* -T), pGC. (11) 

Let us now turn to a rather rare case of the convex FL. Here the first 
intersection between the FL and the CL occurs at the zone diagonals, and 
if one rotates all the zone picture by the angle = 7r/4, these points will 
coincide with the intersection points between the FL and the CL in the case 
of the concave FL. The analysis shows that this feature seems to hold in 
deahng with all the solutions, including the gap landscape. We shall revisit 
this prediction of our model in a future paper. 

So far we have neglected all the electron-electron interactions but the 
longe-range one given by Eq.(^. By involving an electron-phonon exchange, 
the most pronounced out of the remaining ones, we trigger, on one hand, a 
regular component Ar(p) 7^ anywhere in the momentum space. As a result, 
the Landau criterion for superconductivity, violated in the above model at 
T > Tj, is now satisfied. Presumably, the magnitude of Ar(p) slowly varies 
with X, allowing us to estimate it from highly overdoped compounds. Since 
in this case, Tc is small, we infer that the impact of A^ on properties of 
the superfiuid state is insignificant. On the other hand, the electron-phonon 
exchange, specified by the phonon propagator D{uj,k) = A;^/(a;^ — c^/c^), 
gives rise to a renormalization of the constant /„, as well. The respective 
contribution 5A(p) to the gap A is given by the integral 

5A(p) ~ / Dip-p„co,) ,^-|^P;^ dnp: . (12) 

Employing in the integral (|T2|) the "tower" structure of the function Aa(p, x ~ 
Xm), one can decompose ovewhelming contributions into two: one from the 
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same "tower" and the other, from the neighboring one. The first contri- 
bution, proportional to the "tower" range, is smalL When calculating the 
second one, where the momentum transfer p — pi is comparable to pp, the 
propagator D can be replaced by -1, yielding a number, which suppresses the 
initial fa value. This interference may be significant. 

Now we proceed to evaluation of the superfluid density PsiT), expressed in 
terms of a correlator of the velocities 9^°(p) / dp. Evaluation of this correlator 
in crystals with the help of the Larkin-Migdal method [121 El yields: 



Ps{T) - 



2 J dpi 



9n(p,T) dfiE)dap) 



dpi BE dpi 



dr , (13) 



the function n{p, T) being given by Eq.®. In ideal homogeneous Fermi gas, 
where ^(p) = p^ /2M — p, Eq. lfT^ is converted into the ordinary textbook 
formula. Obviously, no contributions to ps are made from regions in momen- 
tum space, where the ratio A(p)/T is negligible. Indeed, if A(p) = 0, the 
distribution n(p,T) is converted to np{p,T) = (1 -t- exp(^(p)/T))^-'^, while 
the product {df{E)/dE){d^{p)/dpi) dnp{p,T)/dpi, and both the terms 
in Eq. (fT^ cancel each other. As a result, at T > A^, contributions from 
regions, other than the C domain, may be neglected. 

A full examination of the formula (fT^ will be done elsewhere. Here 
we restrict ourselves to the region of the dopings x ~ Xm and tempera- 
tures T ^ T* fal'^- In this case, the ratio E{p G C)/T is small, and 
expansion of the terms in the integrand of Eq. lfT^ yields dn{p,T)/dpi ~ 
-(ae(p)/9p,)/(4T) + eip)mp)/^P^)/ilQT^) + A'ip){dap)/dp.)/iA8T') 
and df {E)/dE ~ -l/(4r) + E^ (p) / (IQT^) . After cancelling similar terms 
and employing relation we are left with 

p,(r -> T*) 1 / ^A2(p)^dr ^ an(r - r)V(T*)2 , (14) 
48 Jc opi opi 

where the numerical factor a is of the order of 10^^. As we shall see, such a 
suppression of ps{T — > T*), results in a marked distinction between the criti- 
cal temperature T^. for termination of superconductivity and the temperature 
T* for vanishing of the gap A. The reason for that is a great diversity in the 
gap values, which, as we have seen, results in the temperature shrinkage of 
the domain of integration over momentum space in Eq. ()13|) . 

Strictly speaking, in two-dimensional systems, the temperatures T* and 
Tc never coincide because of the Berezinskii-Kosterlitz-Thouless (BKT) phase 
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transition jlHl UHl 1201 , terminating superconductivity due to a spontaneous 
generation of an infinite number of vortices. This transition always occurs 
before vanishing of the gap A. The BKT temperature, being, in fact, the 
critical temperature Tc, is given by the equation ^Hl 

vrTe = 2p,(T,) . (15) 

In conventional superconductors, where ps{T T*) ~ n{T* — T)/Tc, and, 
hence, (T* — Tc) ~ T^/e'p, the ratio Tc/ep does not exceed 0.2%. However, in 
two-dimensional electron compounds with the doping x ~ x^, the situation 
is different. Indeed, upon inserting the result (jl4j) into Eq. dlSp . one obtains 

T, = a^eUT*-ny{T*y , (16) 

where the factor — 10~^. In high-Tc superconductors, the ratio Tc/e^ is 
of order of 10~^, and, hence, the ratio r = (T* — T^/T* may attain values 
comparable to unity. 

In conclusion, we have evaluated the effect of critical antiferromagnetic 
fluctuations on electron spectra and superfiuid densities of superfiuid states of 
overdoped and optimally doped high-T^ compounds. In underdoped electron 
systems, the situation is more complicated due to the emergence of a branch 
of low-lying collective excitations, whose contribution to properties is signif- 
icant [21] ■ The interplay between these oscillations and critical fluctuations 
in underdoped compounds will be studied a separate paper. 
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